Scaling properties of one-dimensional off-diagonal disorder 
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Validity of the single parameter scaling (SPS) in one dimensional Anderson model with purely 
off-diagonal disorder is being studied. It is shown that the localized region with standard symmetry 
is divided into two regimes: SPS and non-SPS. Scaling relations of the Lyapunov Exponent are 
proposed for these two regimes. In the non-SPS regime, in additional to the localization length, there 
exists a new length scale which is related to the integrated density of states. A physical interpretation 
of the new length is the cross-over length which separates regions with chiral symmetry from those 
that have standard symmetry. 
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I. INTRODUCTION 

It is well-known for about forty years that all electron 
states in standard one dimensional (ID) disordered mod- 
els are localized for any strength of disorder, and there is 
no localization transition in ID systems 1]. However, in 
the case of off-diagonal disorder, there is an anomalous 
localized state at the band center :2iiS i4i5ii . It has been 
proposed that the Lyapunov Exponent (L.E.) is the ap- 
propriate scaling variable to describe fluctuations of the 
conductivity. 

^(N) ^ + = ^-ln(T) (1) 

where g{= T/R) and T are conductance and transmission 
coefficients through the system with length N. 

A reason for the revival of the interest in ID disor- 
dered model is due to the revision of the well-known Sin- 
gle Parameter Scaling (SPS) hypothesis. According to 
this hypothesis there exists a single parameter, con- 
ductance g, which determines scaling properties of g{N). 
Soon after the report of SPS, it became clear that one 
should consider scaling of the full probability distribution 
of conductivity. 

In order to take fluctuations of conductance into ac- 
count, one should consider a parameter 7 (L.E.) defined 
in Eq.(Q instead of conductance g itself. In the thermo- 
dynamic limit, 7o = "f{N — > 00) has a non-random value 
which is the inverse of the localization length A. This 
parameter has normal distribution for iV A, and its 
dispersion a-y obeys the law of large numbers and is not 
an independent variable. 
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where the parameter t is usually defined as a conven- 
tional scaling parameter in literature. The above equa- 
tion was originally derived by Anderson et al. 7] by using 
the random phase hypothesis. 



However, as shown in without the assumption of 
phase randomization, ID SPS (Eq.|2]) is violated, where 
states are much far apart from each other than the local- 
ization length (A). 

This new characteristic length scale £s which is related 
to the distance between states, for the states near the 
band center is defined in terms of total number of states 
{N{E)) whose energy is less than E j^. This crite- 
rion was initially extracted from the exact calculation of 
the variance of L.E. for the Anderson model with Cauchy 
distribution of the site energies M. H owever in contradic- 
tion with violation of SPS in ID Tl], it has been recently 
shown using exact diagonalization 12] and transfer ma- 
trix method that at special point E = 0, SPS holds 
perfectly in ID. 

In 2D case, the interest is motivated by the experimen- 
tal observations of a metal-insulator transition which is 
at odds with the SPS for noninteracting electrons [l^ . 
The validity of SPS in 2D is currently very controver- 
sial. There exists some numerical analysis of 2D Ander - 
son model which confirms the SPS hypothesis [I5ll6|l7|. 
Other studies suggest a two-parameter scaling Il2ll3ll8l | . 
It has been shown 18] that 2D SPS does not follow 
Eq.p. 

In ID systems with off-diagonal disorder (random hop- 
ping model), it is clear that an anomalously localized 
state a.t E = 0, results in a violation of SPS. Divergence 
of the localization length and density of states at the 
band center in this model CHI, is in contradiction with 
the scaling theory. Unus ual pro perties of this model are 
due to chiral symmetry |2ll22j] . In an interval close to 
an anomalous state, SPS does not hold The main 

objective of the present paper, is to answer how far from 
the anomalous state (at E = 0), ID SPS will again be 
held. 

In this paper, with care of some debates on the va- 
lidity of ID SPS, we reexamine the scaling properties of 
one-dimensional system with purely off-diagonal disor- 
der by using transfer matrix method. Our attention is 
on a region near the band center which contains strongly 
localized states with standard symmetry. In the strong 
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localization limit, it will be shown that the L.E. distribu- 
tion function is normal. In this region, it is shown that 
there exists a new length scale which is the same as the 
length scale (Is) defined in Eq.^TTl I- The SPS exists as 
long as the localization length A exceeds is- In the SPS 
region, it is shown that the L.E. only depends on the 
disorder strength as 70 cx o'fn(t)' "^^^ scaling properties 
of the non-SPS region which has been reported in Ref. 
[2^ . is confirmed by a data collapse. The variance and 
mean of the L.E. for different disorder strengths, system 
sizes and also for various range of energy spectrum, lie 
on a single curve when they are expressed in terms of the 
scaling parameter (r) defined in Eq.(|2Jl as a function of 
the ratio k = X/£s- It can also provide a physical inter- 
pretation for £s as a cross-over length between two chiral 
and standard symmetries. 

This article is organized as follows: Section II describes 
our model and the exact calculation of all L.E. moments 
at the band center. Section III describes the transition 
from chiral symmetry region to localized region by the 
calculation of L.E. distribution function and its mean. 
In this Section, it will be shown that the localized region 
is divided into non-SPS and SPS regimes. A new length 
scale(£,j) which controls the scaling theory is defined in 
Section (IV). We try to find a meaningful physical inter- 
pretation for the new length scale as a cross-over length 
in Section (V). Discussions and conclusions are finally 
presented in Section VI. 



II. MODEL AND MOMENTS OF LYAPUNOV 
EXPONENT 

We consider non-interacting electrons in ID disordered 
systems within a tight binding approximation. The 
Schroedinger equation with the assumption of nearest- 
neighbor hopping becomes 



+ U,i+lipi+l + U-l,ii^i 



(3) 



where E is the energy corresponding to the electron wave 
function. is the probability of finding the electron at 
site i, £i are the site potentials and U-i^i — tij^i — U the 
hopping terms. Using the transfer matrix method, one 
can relate the electron wave functions at the two ends of 
the system to each other. In our model, we consider all 
site energies to be zero and a periodic boundary condi- 
tion on hopping terms as ti — tAr+i. All energies which 
appear, are scaled by typical mean of hoppings terms Iq, 
where ln(<o) =< In(ti) >c.a.- Here, c.a. refers to the 
configurational average. The L.E. can be extracted from 
the eigenvalues of the total transfer matrix [2^ . 

As proved in Ref. [2^ . the L.E. at £^ = has a semi- 
Gaussian distribution with a mean which can be derived 
in terms of the pair correlation function. By having the 
distribution function, higher powers of the L.E. can be 
simply derived in the case of correlated and uncorrelated 
disorder at the band center {E = 0) as: 



< 7" - < 7 >^ < 7-^ TT < 7 >■'; 

< 7" >oc< 7 >" (4) 

Therefore, higher moments of the L.E. can be written as: 

=< (7- < 7 >f (J - 1) < 7 >' 



<(7-<7>)' >= (2--)<7>=^ 



(5) 



This can be generalized to the n'th moment of the L.E. 
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FIG. 1: The size dependency of the Lyapunov Exponent for 
different energies with 77 = 0.18. 

which will be as < (7— < 7 >)" >cx< 7 >". By pay- 
ing our attention to the L.E. form at £' = (< 7 >oc 
<^in(t)/N^/'^), it can be seen that the variance of the L.E., 
scales according to the law of large numbers for uncorre- 
lated disorder. As it has been mentioned by Anderson, 
et al 7] , the localization properties can be described by a 
variable (such as L.E.) whose width of distribution func- 
tion follows the law of large numbers. 
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This equation and its equivalence in Eq.jS] are consis- 
tent with the result of Ref. , where it was derived for 
a weak disorder by solving the Fokker-Planck equation. 
However, the size dependence of the L.E. variance will 
change when the disorder becomes correlated. As a re- 
sult, the L.E. distribution function and all its higher mo- 
ments converge for large system sizes. So, L.E. is a good 
variable to describe statistical properties of disordered 
systems. 
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III. SCALING AND DISTRIBUTION 
FUNCTION OF LYAPUNOV EXPONENT 

A. Scaling of Lyapunov Exponent 

We calculate the mean and variance of L.E. by using 
the transfer matrix method when randomness is imposed 
on ln(t)'s. The study of L.E. close to the band center re- 
sults in the coexistence of two symmetries in this system. 
It can be shown that there is a chiral symmetry at i? = 0. 
This is a significant property of purely off-diagonal dis- 
order with nearest-neighbor approximation. However, at 
energies close to the band center, and for lengths greater 
than a cross-over length (Ncr), chiral symmetry is bro- 
ken. At sufficiently long lengths, localization properties 
will flow to those of the standard symmetry class. All 
states in this regime are strongly localized. 

For any realization of the disorder, the energy density 
of states is symmetric around the band center. This sym- 
metry, which originates from the fact that the disorder 
preserves the bipartite structure of the lattice, is referred 
to as chiral symmetry. The chiral symmetry is broken by, 
e.g., on-site randomness or next- nearest- neighbor hop- 
pingim. 

In the case of purely onsite disorder case which was 
originally considered by Anderson one distinguishes 
three universality classes, corresponding to the presence 
or absence of time reversal and spin-rotation symmetry. 
These three classes are called orthogonal, unitary, and 
symplectic |23. Here, we will refer to these as the three 
standard universality classes. In this paper, it will be 
shown that the standard symmetry region (strongly lo- 
calized region) is also divided into two regimes (non-SPS 
and SPS regimes defined in section (III.C) and Fig.Q), 
which depend on the number of scaling parameters. 

Fig.© shows the scaling properties of L.E. near the 
band center. All data with various energies lie on a single 
curve when (7 x E~^) is plotted in terms of the dimen- 
sionless variable {N /Ncr)- It confirms a power law diver- 
gence of the localization length where energies belong to 
the non-SPS regime (Fig.Q). Therefore, the following 
scaling law of L.E. &t E ^ can be proposed p3 |. 



(7 X E-'^) cx 



'^/n(t)(iv;;7) 



iy_N-i/2 



N/N„ « 1 



'ln{t) 



N/Ncr > 1 



where the cross-over length is as 
« E-^'^ /ct 



ln{t) 



(7) 



(8) 



and 77 w 0.18 ±0.03. However, for energies very close 
to the band center, rj has a small energy dependence. 
Fig-O shows a transition from the region with chiral 
symmetry {N <^ Ncr) to the region with standard sym- 
metry {N ^ Ncr)- It was also checked that each of the 
data sets do not collapse on each other when one uses the 
logarithmic energydependence of the localization length 
as seen in Refs.l^la). 



B. Distribution Function of Lyapunov Exponent 

The localization properties of different symmetry re- 
gions can be also characterized by the distribution func- 
tion of L.E. As it was mentioned in section (II), at the 
band center, the distribution of L.E. is semi-Gaussian. 
For zero energy (£■ = 0), numerical evidence in Fig.(|21a) 
confirms such analytical distribution function for all sys- 
tem sizes. 
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FIG. 2: Lyapunov Exponent Distribution Function for differ- 
ent system sizes at a) i5 = and h) E = 10~^. Disorder 
strength is considered to be 0.1. The number of samples is 
2 X 10^ configurations. 



For energies near the band center {E ^ 0), and for sys- 
tem sizes smaller than the cross-over length (A^ ^ Ncr), 
the distribution of L.E. is semi- Gaussian. In this region, 
the band center behavior is dominant. However, as the 
system size increases, the distribution function becomes 
more Gaussian-like for N ^ Ncr- This region has a 
regular Anderson-like behavior with standard symmetry 
class. Fig.(|33) shows distribution function of the L.E. 
for the same system sizes as the band center case. It 
can be seen that there exists again a transition between 
these two symmetries; from chiral to standard symme- 
try or from semi-Gaussian to Gaussian distribution. For 
the sake of completeness, the skewness of the distribu- 
tion functions has been calculated as a measure of the 
symmetry of the distribution. It can be defined as p^: 



Skewness 



< (7 - 7)3 > 
< (7 - 7)^ > 



3/2 



(9) 



Distribution close to the normal form, has a skewness 
equal to zero. A distribution whose skewness has ab- 
solute value less than 0.5 is considered fairly symmetri- 
cal. Therefore, distributions of long enough systems in 
Fig.(I2b)i are very close to the Gaussian form. In the 
SPS regime (Fig.©) where the L.E. is independent of 
size and energy, distribution function of L.E. is exactly 
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Gaussian (the skewness order of 10~^) and independent 
of size. All distribution curves have been softened by 
the Kernel smoothing method jl^ without changing any 
statistical characteristic of distributions. 



C. SPS and Non-SPS Regimes 

For a fixed length (N), Eq.(|2j) proposes a critical en- 
ergy point which separates two regions with different 
symmetries. In fact, for energies greater than ei]-'. cx 
(TVcTj^^^j^ ) 2^ , the system is in the localized region. Fig. © 
shows the energy dependence of L.E. in the localized re- 
gion. It can be seen that there is a second critical point ( 

eir^ w 10~^ in Fig.Q) where the energy spectrum is di- 
vided into SPS and non-SPS regimes. In the SPS regime, 
L.E. is independent of energy and the scaling theory is 
vahd (Eq.Q). Although, near the band edges anomaly, 
L.E. will become energy dependent. In the SPS regime, 
it has been checked (Figs.lPEl)) that in contradiction to 
the result of Rcf. 27], the L.E. only depends on disorder 
strength. Fig.Q shows that the L.E. is proportional to 
the square of disorder strength. 



70 OC cri„(() 



(10) 



The line Fitted on data in Fig.l^J has a slope equal to 
the value 2. The coefficient of the above scaling law 
is ( 3 oo±o 2 ) "^hich can be extracted for a fixed disorder 

strength (crin(t) = 0.1) and for an energy (e = 0.1 > £cr.) 
in the SPS regime. The skewness of this point is about 
0.005. This skewness shows an exactly Gaussian form for 
L.E. in this regime. 

It can be seen that matching of two Eas. H7llO() at the 

(2) 

boundary of non-SPS to SPS regime {E — Scr ) leads to 

f 2) 

a second critical point (scr ) which is independent of all 
system parameters (a constant). 

Now, we numerically study the variance of L.E. as a 
fimction of system parameters. Fig© shows Na?^ versus 
energy at fixed system size. As it can be seen, the vari- 
ance of L.E. is approximately independent of energy at 
energies in the SPS and also non-SPS regimes. Figs. (|3l4() 
show that size and disorder strength dependence of the 
variance of L.E. follows from Eq.®. As it is clear from a 
fitted line (with slope 2) on data in Fig.Q), the quantity 
Na?^ is proportional to fi^n(t)- ^^^^ dependence of 

the L.E. variance (CT^(iV)) is shown in the inset Fig.Q). 
The L.E. variance decreases with the inverse of the sys- 
tem size similar to the size dependence of variance at 
the band center (Eq.®). So, the variance of L.E. at the 
band center can be generalized to other energies near the 
band center. 



IV. VIOLATION OF SINGLE PARAMETER 
SCALING 

According to Eq. (O , the two parameters of the distri- 
bution reduce to only one. Parameter r can be defined as 
a measure of SPS in that equation. Now, in this Section, 
we try to find the independent parameters of the system 
in non-SPS regime. It is again stressed that the random 
hopping model away from the E = and in the localized 
region is being studied. The quantity of interest is £s 
that is related to the integral density of states from the 
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FIG. 3: SPS to non-SPS transition with plotting Lyapunov 
Exponent and Ncr?^ versus energy for system size 10^. eil^. 
which depends on the disorder strength, is about 5 x 10~^ 
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FIG. 4: Standard deviation of the Lyapunov Exponent versus 
disorder strength for system with length 10^ . Dependency of 
70 to disorder strength in the SPS regime {E = 0.1). Inset 
figure shows the system size dependence of the variance of 
Lyapunov Exponent. 
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band center to a given energy, normalized by the total 
number of states in the band. 



(11) 



sm{nN{E)) 



For the Anderson model, N(E) can be computed by 
the node-counting theorem [29|. By starting an initial 
vector in transfer matrix method, we count the number 
of wave function nodes as the length is scanned. 

Fig-© shows the numerical result of inverse scaling 
parameter (t(k)) in terms of the dimensionless parameter 
(k = J-) for different values of disorder strength and 
energy. 

The data included in this graph correspond to the lo- 
calized regime with standard symmetry class where the 
L.E. has a Gaussian distribution. The inset Fig.jSJ) shows 
that the skewness of the data are less than 0.25 for 
K > 0.1. All data are in the region N ^ £s. What is 
important, is that all data with different values of energy 
and disorder strength and also system size collapse to a 
single curve when they are expressed in terms of 1 /r and 
K. Therefore, for k <C 1 (non-SPS regime) variance of 
L.E. depends on two parameters; k and the mean of L.E. 
In the case of k ^ 1 , the inverse of the scaling parameter 
(l/r) in the present model goes to the value 0.5 which is 
different from unity. Therefore, our expression from non- 
SPS is only the deviation of Eq.Q from 0.5. Energies are 
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be seen that for k ^ 1, SPS is clearly violated (non-SPS 
region). For k 3> 1, independent of disorder strength, 
standard SPS is restored again. In the SPS spectral re- 
gion, the localization length is a macroscopic length as 
only one parameter of system. The conductance can be 
defined as a function of this variable. 

Both length scales is and the localization length are 
decreased when energy is swept from the vicinity of the 
band center to the band edges. In energies close to 
zero, the new length scale is greater than the localiza- 
tion length. There is a critical energy where both length 
scales are of the same order (A « Ig). Far from the 
band center, the localization length is independent of en- 
ergy and the new length scale steeply decreases so that it 
would be much smaller than a macroscopic localization 
length. 

This result confirms the general conjecture of the au- 
thors that the second moment of the distribution func- 
tion of L.E. can be universally described in terms of vari- 
ables T and K regardless of the microscopic nature of 
the models under consideration. The form of the func- 
tion t(k) may differ for different models and its essential 
behavior is not universal. All models follow t{k) = 1 
for K ^ 1, while in the hopping disorder model, it is 
t(k) ~ 2. In the model studied in the present paper, 
for K <^ 1, it can be seen an exceptional behavior com- 
pared to other models such as Lloyd model .8„ ,10., . 2_3] and 
Anderson (onsite disorder) and superlattice models [23l| . 
The scaling parameter t increases with k in the hopping 
disorder model for k <^ 1, while in the above models, 
T steeply decreases with k. As an example, analytical 
calculations carried out in Ref. ^ for the Lloyd Model 
produced r — in /2)k. 

The power law form is the best fitted curve for k ^ 1. 
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(12) 



Coefficients are estimated by using a linear regression in 
log-log plot. The fitted power of k and its coefficient are 
P = 0.815±0.008 and a = 1.64±0.03, respectively. Since 
there is a kind of delocalization at the band center, the 
L.E. sharply decreases near the band center compared to 
Na?^ (Fig.©). Therefore, the inverse scaling parameter 
decreases in energies close to the band center. 



FIG. 5: Measure of the validity of the Single Parameter Scal- 
ing by inverse Scaling Parameter (l/r = jo/{Na?^)) in terms 
of K = X/£s for disorder strengths from crin(t) = 0.05 to 0.2. 
Energies depending on disorder strength, are scanned from 
E = 10"'^ - 10°. Here N '> is. Inset figure shows the skew- 
ness of the data. The number of samples is 10* configurations. 

scanned from near the band center to the middle of the 
conduction band, and far from the second anomaly in the 
band edges. In fact, close to any anomaly (here a delocal- 
iztion at the band center), a violation of SPS will happen 
llOl in an interval close to the anomalous state. It can 



V. CROSS-OVER LENGTH AS A PHYSICAL 
INTERPRETATION OF is 

Numerical results in Fig.||2Jl show that as a crude ap- 
proximation, one can consider the variance of L.E. inde- 
pendent of energy in the non-SPS and SPS regimes (not 
near their transition point). It was shown that the form 
of variance near the band center is similar to its form 
in the band center (Eq.®). By using this form of the 
variance and the scaling relations derived in Eas. (|7ll0|l . 
the scaling parameter function (r) can be proposed in 
the non-SPS and SPS regimes. 
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First, we investigate the scaling parameter in the non- 
SPS regime. According to the infinite L.E. in the local- 
ized region (Eq.Q) 70 oc (J^^^^^E^ and its variance as 
it can be seen that the measure of SPS 



Nal oc cr 



.2 

'7 ln(t)' 

(scaling parameter r) can be proposed to have a linear 
relation with the dimensionless variable k' — X/N^r- 



- oc i;" oc — 
T 7V„ 



(13) 



This equation shows a deviation from the SPS value 
(unity) in the non-SPS regime. The cross-over length 
Ncr plays the role of a length scale like £s in this system. 
The above form for the scaling parameter is independent 
of the system parameters such as disorder strength when 
it is expressed in terms of r and k'. This expression 
is confirmed by Fig.J^)) which shows the inverse scaling 
parameter versus k'. It can be seen that for k' ^ 1, all 
data for different disorder strengths, coincide with each 
other on a single curve. However, since in the second 
critical point (k' « 1), the variance of L.E. is energy 
dependent and also, the scaling form of L.E. (Eq.©) is 
not correct in this point (Fig. (PI), curves with different 
disorder strengths are separated from each other in the 
transition point. 

In the SPS regime, the scaling of L.E. proposed in 
Ea. (|10|l 7o oc cj^nj-t) and its variance form as Eq.®, show 
that the scaling parameter r is independent of system 
parameters and a constant (Fig.lO). 

Therefore, both the new length scale and the cross-over 
length that are of the same order in the non-SPS regime, 
can characterize the scaling properties. Fig.Q compares 
these two length scales. It shows the ratio of the new 



length to the cross-over length in terms of energy. In 
the Non-SPS regime, they both weakly depend on en- 
ergy, although, in the transition region, the difference is 
remarkable. 

As a result, the cross-over length is proposed as a phys- 
ical and meaningful interpretation for the new length 
scale (^s)- On the other hand, we showed that the sta- 
tistical distribution of L.E. is different for sizes lower or 
greater than this length scale. 
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FIG. 7: The ratio of the new length scale {£3) to the cross- 
over length in terms of energy for different values of disorder 
strength. 



VI. CONCLUSION 
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FIG. 6: Inverse scaling parameter 1/r = 'yo/{Na'^) versus 
localization length in the ratio of the cross-over length {k = 
^/Ncr) for different disorder strengths. Energies depending 
on disorder strength, are scanned from E = 10~^ — 10°. 



In this paper, we study the scaling theory in the hop- 
ping disorder model. The main result of this paper is 
to show the single parameter scaling (SPS) is violated 
not only in a region with chiral symmetry, but also in 
the localized region where there exists a standard sym- 
metry class. The localized region is also divided into two 
regimes: SPS and non-SPS regimes. We proposed the 
scaling relations for the Lypunov Exponent in these two 
regimes. The criterion of the SPS is controlled by a new 
length scale which is related to the integral of density of 
states, is defined in Ref.Q. The SPS holds when the 
localization length A exceeds the new length (A ^ £s)- 
In A <^ £s regime, standard deviation of the Lyapunov 
Exponent ( 7) distribution can be described by two inde- 
pendent scaling parameters: the mean of 7 and k = X/ Is- 

We showed that all data related to the variance and 
mean of the Lyapunov Exponent with different values 
of disorder strengths, system sizes and also the data ex- 
tracted from various energy regions, lie on a single curve, 
when they are expressed in terms of the inverse scaling 
parameter 1/r = 70 /Na"^ and the dimensionless variable 

At. 
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The cross-over length (Ncr) which separates the region 
with chiral symmetry from that of standard symmetry, is 
proposed as a meaningful physical interpretation for £s- 
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